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Abstract 

In this paper we have proposed an almost unbiased estimator using known value of some 
population parameter(s). Various existing estimators are shown particular members of the 
proposed estimator. Under simple random sampling without replacement (SRSWOR) scheme the 
expressions for bias and mean square error (MSE) are derived. The study is extended to the two 
phase sampling. Empirical study is carried out to demonstrate the superiority of the proposed 
estimator. 

Key words: Auxiliary information, bias, mean square error, unbiased estimator, 
two phase sampling. 

1. Introduction 

Consider a finite population U = Ui, U 2 , ...Un of N units. Let y and x stand for the 
variable under study and auxiliary variable respectively. Let (yj, Xj), i=1, 2,.., n 
denote the values of the units included in a sample s n of size n drawn by simple 
random sampling without replacement (SRSWOR). The auxiliary information has 
been used in improving the precision of the estimate of a parameter (See 
Cochran (1977), Sukhatme et. al. (1984) and the references cited there in). Out 
of many methods, ratio and product methods of estimation are good illustrations 
in this context. 

In order to have a survey estimate of the population mean Y of the study 
character y, assuming the knowledge of the population mean X of the auxiliary 
character x, the well-known ratio estimator is 




Bahl and Tuteja (1991) suggested an exponential ratio type estimator as - 



(1.1) 
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Several authors have used prior value of certain population parameter(s) to find 
more precise estimates. Sisodiya and Dwivedi (1981), Sen (1978) and 
Upadhyaya and Singh (1984) used the known coefficient of variation (CV) of the 
auxiliary character for estimating population mean of a study character in ratio 
method of estimation. The use of prior value of coefficient of kurtosis in 
estimating the population variance of study character y was first made by Singh 
et. al. (1973). Later used by Singh and Kakaran (1993) in the estimation of 
population mean of study character. Singh and Tailor (2003) proposed a modified 
ratio estimator by using the known value of correlation coefficient. Kadilar and 
Cingi (2006), Khosnevisan et. al. (2007), Singh et. al. (2007) Singh and Kumar 
(2009) and Singh et. al. (2009) have suggested modified ratio estimators by 
using different pairs of known value of population parameter(s). 

In this paper under SRSWOR, we have proposed almost unbiased estimator for 
estimating Y. 

2. Almost unbiased ratio type estimator 



Such that t 0 t rB ,t rEe €w r where w r denotes the set of all possible ratio type 
estimators for estimating the population mean Y. By definition the set w r is a 
linear variety, if 



where ui i (i=0, 1 , 2) denotes the statistical constants and R denotes the set of 
real numbers. 

To obtain the bias and MSE of t w , we write 

y = Y(l + e 0 ), X = x(l + e 1 ), 

such that 

E (e 0 ) =E (ei)=0. 

E(e§) = fjC* , E(e?) = i x c \ , E(e oei ) = f lP C y C x . 



Suppose 

t0 




(aX+b)“(ax+ b 
(a_X-i-b) + (ax+ b 



i) 



tvvr = “<>y + «l*r. + “ztrse ^ 
for =1, <*>;eR 



(2.1) 

(2.2) 
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Expressing t w in terms of e’s, we have 

t w = Y(1 + e 0 )[a> 0 + 0^(1 + ee^ -1 + trt z exp {- ^{1 + Oe^ -1 }] 



(2.3) 



where 0 = 



aX 

aX+b' 



Expanding the right hand side of (2.3) and retaining terms up to second order of 
e’s, we have 



t — Y 

L w — 1 



1 4“ 6q — to0e-j "f 0 2 "I - ej — Bue^e^j 



where = (o^ + =). 



(2.4) 

(2.5) 



Taking expectations of both side of (2.4) and then subtracting Y from both side, 
we get the bias of the estimator t w , up to the first order of approximation as 

B (‘„) = f l 7 (», + Jf) - 9<*)pC y C x ] (2.6) 



B(t„) = f,Y[S 3 C. ! - OpC y cJ 



B(t rs =) = f,Y 



30 3 c£ 



y 

6pCyC X 



From (2.4), we have 

(t w -Y) Y[ eo -0coej 



(2.7) 

(2.8) 

(2.9) 



Squaring both sides of (2.9) and then taking expectations, we get MSE of the 
estimator t w , up to the first order of approximation, as 

MSE(t w ) = fjY [C;. + 0 2 co 2 C| - 20ojpC y Cj (2.1 0) 

This is minimum when 

t o = k(=p^). (2.11) 

Putting this value of u>(= k) in (2.10), we get the minimum MSE of t w as 
min. MSE(tJ = fJ z C z (l- p 2 ) (2.12) 



which is same as that of traditional linear regression estimator 
from (2.5) and (2.11), we have 

^ + ^ = (2.13) 

From (2.2) and (2.13), we have only two equations in three unknowns. It is not 
possible to find the unique values for to [ S ,i = 0,1,2, In order to get unique values 
for to iV we shall impose the linear restriction 



co 0 B(y) +^B(t rB ) +o> 2 B(t rse ) = 0 (2.14) 

Equations (2.2), (2.1 1) and (2.14) can be written as in the matrix form as 



rl 


1 


1 1 


-OV 




T 


0 


l 


1/2 


(Oj 


= 


k 


.0 


B(t rs ) 


B(t rge )J 


“2- 




0- 
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Using (2.1 5) ,we get unique values of o> i B (i=0,1 ,2) as 

“o = £ 

©i = — 

1 A r 
^2 

OJ? = — 

2 Ar' 

where 

A r = B(t rse )-i B(t„) 

A r0 = B(t rse ){l - k} B(t Pi ) {k- 

A r= k-BCO 

A r2 = 



(2.16) 



(2.17) 



Use of these (i=0,1 ,2) remove the bias up to terms of order (n' 1 ) at (2.1 ). 



3. Product -type estimators 

Suppose t 0 = y , t p! =yggj), 



L pse 



_ f (ax ■+■ b‘i - ( aX-H b') 1 

= ^ l (ax+b) + (aX+bJ 



(3.1) 

(3.2) 



such that t 0 , t pg/ t pse e Q, where Q denotes the set of all possible product -type 
estimators for estimating the population mean Y. By definition, the set Q is linear 
variety if 

= q &y + qit pB + q 2 t PB e e Q 

for z? =0 q i = l, qj e R 

where qj(i=0,1 ,2) denotes the statistical constants. 

Expressing t q in terms of e’s, we have 

t q = y(l + e o) [to D + (1 + 0e t ) 4- o) 2 exp (l + 06!) _1 }] 

where 6 = -4^ 



(3.3) 



aX+b 



Expanding the right hand side of (3.3) and retaining terms up to second power of 
e’s, we have 



t q = Y [l + e 0 + 0qo, - ej + gSe^] 



where q = qj + y 



(3.4) 

(3.5) 



Taking expectations of both sides of (3.4) and then subtracting Y from both 
sides, we get the bias of the estimator t q , up to the first order of approximation 
as 



[-fe 2 <2+qe P c y c, 



(3.6) 
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Bias expression for the estimators t pE and t pBe is given by 

bW = t J [epc/J 

From (3.4), we have 
(t,-Y3^Y[e o +0q ei ] 



(3-7) 

(3.8) 

(3.9) 



Squaring both the sides of (3.9) and then taking expectations, we get MSE of the 
estimator t q , up to the first order of approximation, as 

MSE(t q ) = f\Y 2 [cj. + 0 2 q 2 C 2 4- 20qpC y C x ] (3.1 0) 

which is minimum for 

q=-k=-p^ (3.11) 

Putting this value of q(=-k) in (3.10), we get the minimum MSE of t c as 

min. MSE(t q ) = f 1 Y 2 C 2 (l - p z ) (3.12) 

which is same as that of traditional linear regression estimator. 

From (3.5) and (3.1 1 ), we have 

q,+^ = -k (3.13) 



From (3.2) and (3.13), we have only two equations in three unknowns. It is not 
possible to find the unique values for qj’s, i=0,1 ,2. In order to get unique values of 
qi’s, we shall impose the linear restriction 

q 0 B(y) +q 1 E(t ps ) + q 2 E(t pse ) = 0 



pey 1 'iZ^'Y pee 

Equations (3. 2), (3. 13) and (3.14) can be written in the matrix form as 

1 1 1 
0 1 1/2 

0 B (v) B ('p..) 

Solving (3.15), we get the unique values of qi’s (i=0,1 ,2) as- 

•Apo’x 









1 ‘ 




9i 


= 


-k 




- t l2- 




- 0 - 



(3.14) 



(3.15) 



q 0 = 



qi = 



Ap 

Yji 



q, = ^ 
V 

where 



(3.16) 
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A p = B(t pse )-i B(tJ \ 

A p o— B (t pse Xl + k} + B(t ps ) j— k — jj > (3.17) 

A pl= “ k - B ( t p g e) 

A p2 = ^ B ( t pe) > 

Use of these qi’s (i=0,1 ,2) remove the bias up to terms of order o(n _1 ) at (3.1 ). 

In Appendix A we have listed some of the important known estimators of the 
population mean, which can be obtained by suitable choice of constants 

CL> if i = 0 , 1,2 , qj, i = 0 , 1,2 and a and b. 



4. Proposed estimators in two phase sampling 

When X is unknown, it is sometimes estimated from a preliminary large sample of 

size n’ on which only the characteristic x is measured (for details see Singh et. al. 
(2007)). Then a second phase sample of size n (n < n’) is drawn on which both y 
and x characteristics are measured. Let x = denote the sample mean of 

x based on first phase sample of size n’, y = and x = ^ be the 

sample means of y and x respectively based on second phase of size n. 



In two phase sampling the estimator t wr will take the following form 

twd = “MY + “id*™! + ^dtrsed e w d (4-1) 

for Zf=o W id = 1, d) id ER (4.2) 

_/aX+b\ _ f (aX +b|— (ax+bl 

where *« = y LsJ and *»• = y ex p 

To obtain the bias and MSE of t wd , we write 
y = Y(l + e 0 ), x = X(l + e[), x = X(1 + 
such that 

E (e 0 ) = E ( ei ) = E( ei ) =0. 

E(4 ) = f l C 2 y , E(ef ) = fjC^ , E(ef ) = f 2 C 2 x E^ ) = fipC y C x 



E(e 0 e' 1 ) = f 2 pC y C x E(e 1 e' 1 ) = f 2 C 



where 



fi =(--^) 
n N 



f2 ( n’ N } 



Following the procedure mentioned in section 2 and 3, we get bias and MSE of 



twd as 

B (t„d) = y [e 2 C, 2 f 3 (u>,d + ^) - Bp C y C,f 3 (<b ld + ^)] (4.3) 
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MSE(t ud ) = Y 2 [fjCy + — 2f 3 pC C x to d ] 



(4.4) 



where f 3 = - — = f-, — f 2 . 

J n n 1 z 

MSE (t Md ) is minimum, when 






= 03 d = k. 



(4.5) 



Putting this value of o> d in (4.4), we get the minimum MSE of t^ d as 

min. MSE(t^ d ) = Y Z C Z [f, - f 3 p z ] (4.6) 

This is same as that of traditional two phase linear regression estimator. 

The bias expression for the estimators t rsd and t rsde is respectively given by 

= — (4-7) 



B(t rBd6 ) = Y 9 2 C|f 3 



3 



— 0pC y C x f 3 -j 3 - 



(4.8) 



From (4.2) and (4.5), we have only two equations in three unknowns. It is not 
possible to find the unique values for w^’s i=0,1 ,2. 



In order to get unique values of w icl , we shall impose linear restriction 

“od B(y) + ti) ld B(t rBd ) + w zd B(t rsed ) = 0 

Equations (4.2), (4.5) and (4.9) can be written in matrix form, as 



rl 


1 


1 1 


■“od‘ 




T' 


0 

.0 


l 

B(t rsd ) 


1/2 

B(W)J 


“id 

i“2d- 


= 


k 

O- 



(4.9) 



(4.10) 



Solving (4.10), we get the unique values of a> id s r ( i = 0,1,2) as 



lod^ 

■*Vd 
A m 

^rd 

r'i — ^ 2d 
0)211 “ ^ 



“od = 



“id = 



(4.11) 



where 

^rd = B(t Pgsd ) - 1 B(t pg( j} 

^rld — ^^(^rsed) 



^r2d — ■BCtj.jd) 



(4.12) 



Use of these Wid’s (i=0,1,2) will remove the bias up to terms of order 0(n _1 ) at 
(4.1). 
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The estimator t written in (3.1), in two phase sampling, will take following form 

tqd = aw? + qidtpsd + qzdtpsed e Qd (4.1 3) 

For £? =0 q Ld =l, Qid eR - (4.14) 



where c? £<i (i=0,1 ,2) denotes the statistical constants . 
The estimators t psd and t psde are 



*psd Y 




and 



t 



psed. 



_ ( Cax+b)-faX - +b1 

l Cax+b)+(aX'+b ] 



Following the procedure of section 4, we get the unique values of qid’s (i=0,1 ,2) 
as 



Qod = 
qm = 

^2d = 



■^pod. ^ 
-^pd 
^pid ^ 

■ipd 

^p2d 



(4.15) 



where 

Apd- BCtrsed) - J B ( t rsd) 

A p od= B(t rsed ){l + k> +| B(t red ){-k-i] ^ 

4pid = — ^ B (t paed ) 

A p2d = k. B (l: psd ) J 



where 

B(t pgd )=Y[0q ld f 3 pC y Cj 

B (‘ P .,d) = f 3 pc y c A -^q ;d f 3 cj' 



The minimum MSE of t qd is given by 



(4.16) 



(4.17) 

(4.18) 



MSE(t qd ) = Y 2 Cj[q-f 3 p 2 ]. 



5. Empirical study 

For empirical study we use the data sets earlier used by Kadilar and Cingi (2006) 
(population 1) and Khosnevisan et. al. (2007) (population 2) to verify the 
theoretical results. 
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Data statistics 



Population 


N 


n 


Y 


X 




c. 


P 


P2OO 


Population 1 


106 


20 


2212.59 


27421 .7 


5.22 


2.10 


0.86 


34.57 


Population 2 


20 


8 


19.55 


18.8 


0.355 


0.394 


-0.92 


3.06 



Table 5.1 : Values of and <?j’s 



CO 1 9 s 


Population 1 


Population 2 


&>Q 


8.590718 


7.892148 


(<lo) 


(21.417) 


(2.919085) 




11.86615 


5.234461 


<0 


(16.14158) 


( 3.576773) 




-19.4569 


-12.1266 


(<lz) 


( -36.5586) 


( -5.49586) 



The percent relative efficiencies (PRE) of various estimators of Y are computed 
and presented in Table 5.2 below. 



Table 5.2: PRE of different estimators of Y 



Estimator 


PRE (Pop 1) 


Estimator 


PRE (Pop II) 


to 


100 


qo 


100 


ti 


212.816 


qi 


526 


t2 


212.803 


q2 


550.261 


t3 


212.606 


qs 


645.256 


t 4 


212.815 


q4 


534.592 


t5 


212.716 


C|5 


581.732 


t6 


212.810 


qe 


465.501 


t 7 


143.992 


qy 


384.447 


t8 


143.923 


qs 


285.920 


t9 


143.988 


qe 


338.487 


tio 


143.990 


qio 


374.584 


tn 


143.991 


qn 


345.118 


1 1 2 


143.959 


qi2 


231.602 


1 1 3 


143.991 


qi3 


424.194 


tl4 


143.987 


qi4 


360.086 


tl5 


143.992 


qis 


356.520 


tl6 


143.911 


qie 


467.051 


tw(0pt) 


384.025 


tq(opt) 


650.263 
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In order to see the performance of the suggested estimators in two phase 
sampling we use the data set of Murthy (1967) (Population III) and Steel and 
Torrie (1960) (Population IV). 



Population 




Cm 


P 


N 


n’ 


n 


Population 1 


0.3542 


0.9484 


0.9150 


80 


30 


10 


Population2 


0.4803 


0.7493 


-0.4996 


30 


12 


4 



Table 5.3 : The values of <*> id '& and q id 's 



"t/‘ s 

(Qid's) 


Population III 


Population IV 


10 Od 


-0.241523 


3.011435 


(Qod) 


(1.808833) 


(1.089979) 




-0.558071 


1 .370950 


(Q Id) 


(0.125381) 


(0.730464) 


w 2d 


1 .799595 


-3.382385 


(^d) 


(-0.934214) 


(-0.820443 ) 



The percent relative efficiencies of various estimators of Y in two phase sampling 
are computed and presented in Table 5.4 below. 



Table 5.4 : PRE of different estimators of Y in two phase sampling 



Estimator 


PRE (Population 1) 


PRE (Population II) 


y 


100 


100 


^rsd 


36.642 


24.562 


tpsd 


4.849 


59.770 


Vsde 


200.420 


48.365 


^psed 


23.628 


115.142 


tud 


276.156 


63.452 


t qd 


34.321 


123.762 


topt 


276.156 


123.762 



6. Conclusion 

From theoretical discussion and empirical study we conclude that the proposed 
estimators under optimum conditions perform better than other estimators 
considered in the article. The relative efficiencies of various estimators are listed 
in Table 5.2 and 5.4. 
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Appendix A 



Table A.1 : Some members of the proposed family of estimators - 



a 


b 


Cflo" 


ft)l 

Oil) 


(q 2 ) 


Ratio Estimator 
(corresponding to 

Wfj=0,1,2) 


Product Estimator 
(corresponding to q ( ,i=0,1,2) 


0 


0 


1 


0 


0 


t 0 = y 

The mean per unit estimator 


q c = y 

The mean per unit estimator 


1 


0 


0 


1 


0 


h = y~- 

The usual ratio estimator 


q!=y| 

The usual product estimator 


1 


c* 


0 


1 


0 


X+C x 

t 2 =y_ * 

X +C x 

Sisodia and Dwivedi (1981) 
estimator 


x+C x 

q 2 =y— 

X+C x 

Pandey and Dubey (1988) 
estimator 


1 


P2OO 


0 


1 


0 


_X + P 2 (x) 

= V 

x + p 2 (x) 

Singh et. al. (2004) estimator 


_ x + p 2 (x) 

= y= 

J X + P 2 (x) 

Singh et. al. (2004) estimator 


p 2 oo 




0 


1 


0 


Xp 2 (x) + C x 

14 = y 

xp 2 (x) + C x 

Upadhyaya and Singh (1999) 
estimator 


-Xp 2 (x) + C x 
q 4 = y — 

y Xp 2 (x) + C x 

Upadhyaya and Singh (1999) 
estimator 




P 2 00 


0 


1 


0 


t _.XC x+ p 2 (x) 

5 xC x +p 2 (x) 

Upadhyaya and Singh (1999) 
estimator 


- xC x +P 2 (x) 

q 5 = y^^ 

y XC x + P 2 (x) 

Upadhyaya and Singh (1999) 
estimator 


1 


p 


0 


1 


0 


rx + pi 

|x + pj 

Singh and Tailor (2003) 
estimator 


_-F + Pl 

qs 7 li + pJ 

Singh and Tailor (2003) 
estimator 


1 


0 


0 


0 


1 


- rx-xi 


- rx-xi 

rt - , — \r pyty 


1 7 - y ex P v , - 

[X + x 

Bahl and Tuteja (1991) 
estimator 


hi - y ex P — _ 
[X + x J 

Bahl and Tuteja (1991) 
estimator 
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0>- 



(q 0 ) 



(qj 



(q 2 ) 



Ratio Estimator 
(corresponding to 
j=0,1 ,2) 



Product Estimator 
(corresponding to q,-j=0,1,2) 



P2OO 



r 8 =yexp 



X-x 



X+x + 2/? 2 (x) 
Singh et. al. (2007) estimator 



x-X 

% = jex^ 

X +X +2yfi, (a) 
Singh et. al. (2007) estimator 



t 9 = yexp 



X-x 



X + x + 2 C, 



Singh et.al. (2007) estimator 



x-X 

<i 9 =yexp 

X + x + 2C ; 
Singh et.al. (2007) estimator 



X-x 






x-X 


X + x + 2 p 




qio - yexp 


X + x + 2 p 



1 10 =yexp 



Singh et. al.(2007) estimator 



Singh et. al.(2007) estimator 



P2O) 



/3,(x)(X-x) 



t ,=yexp 

\fr(x)(X+x)+?C x _ 

Singh et. al. (2007) estimator 



q t =yexp 



fi(x)(x-X) 



/i(x) QC+x)+'2C x 
Singh et. al. (2007) estimator 



P2OO 



L Q(X-x) 

1 2 =^exp — 

\C x (X+x)+2/i(x) 

Singh et. al. (2007) estimator 



^ C x (x-X) 

^12 



C x (X+x)+2&(x) 
Singh et. al. (2007) estimator 



C x (X—x) 






C,(x-X) 


C x (X -\-x) + 2 p 




^13 — y ex p 


C x (X + x) + 2 p 



?i3=yexp 



Singh et. al. (2007) estimator 



Singh et. al. (2007) estimator 



fXX-x) 






p(x-X) 


fXX -\-x) +2C7. 




<?14 “ J eX P 


p{X + x) + 2 C x 



?i4=yexp 



Singh et. al. (2007) estimator 



Singh et. al. (2007) estimator 



P2OO 



?i 5 =yexp 



(X(x)(X-x) 



jX (a) (X +A') +2yfA 

Singh et. al. (2007) estimator 



g 15 =yexp 



J%(x)(x-X) 



/i(x)(X+x)l 2/7 
Singh et. al. (2007) estimator 



P2OO 



f(X—x) 



?i^=3 ; exp — — 

_p(X+x)+2yQ,(x) 
Singh et. al. (2007) estimator 



fix-X) 

<7i 6 = yexr = 

_p(X4X)+2/3,(x) 

Singh et. al. (2007) estimator 



74 



Pak.j.stat.oper.res. Vol.IV No.2 2008 pp63-76 



Almost Unbiased Estimator for Estimating Population Mean using Known Value of Some Population Parameter(s) 



In addition to above estimators a large number of estimators can also be 
generated from the proposed estimators just by putting different values of 
constants io L ,i = 0,1,2 , q t ,i = 0,1,2 and a and b. 
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